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Abstract 

In this paper the general plane strain problem of adhesively bonded 
structures which consist of two different orthotropic adherends is con- 
sidered. Assuming that the thicknesses of the adherends are constant 
and are small in relation to the lateral dimensions of the bonded region, 
the adherends are treated as plates. Also, assuming that the thickness 
of the adhesive is small compared to that of the adherends, the thickness 
variation of the stresses in the adhesive layer is neglected. However, 
the transverse shear effects in the adherends and the in-plane normal 
strain in the adhesive are taken into account. The problem is reduced 
to a system of differential equations for the adhesive stresses which 
is solved in closed form. A single lap joint and a stiffened plate under 
various loading conditions are considered as examples. To verify the 
basic trend of the solutions obtained from the plate theory and to give 
some idea about the validity of the plate assumption itself, a sample 
problem is solved by using the finite element method and by treating the 
adherends and the adhesive as elastic continua. It is found that the 
plate theory used in the analysis not only predicts the correct trend for 
the adhesive stresses but also gives rather surprisingly accurate results. 
The solution is obtained by assuming linear stress-strain relations for 
the adhesive. In the Appendix the problem is formulated by using a non- 
linear material for the adhesive and by following two different approaches. 


1. Introduction 

Generally an adhesively bonded structure consists of three compo- 
nents of different mechanical properties, namely the two adherends and 
the adhesive layer. Even though under most operating loads and environ- 
mental conditions the adherends may behave in a linearly elastic manner. 


(*) This work was supported by NASA-Langley Research Center under the 
Grant NGR 39-007-011 and by NSF under the Grant ENG 78-09737, 
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under relatively severe loading and temperature the adhesive may exhibit 
viscoelastic and/or nonlinear behavior. However, because of the material * 
nonhomogeneity and the geometric complexity of the medium, the exact 
analytical treatment of the related structural problem is hopelessly 
complicated. The existing analytical solutions have, therefore, been 
carried out under certain simplifying assumptions in formulating the 
problem. The primary factors influencing the choice of a particular 
idealized model for the adhesive and the adherends appear to be the 
adhesive-to-adherend and adherend-to-adherend thickness ratios and the 
ratio of the adherend thickness to the lateral joint dimensions. Thus, 
in [1,2] the adhesive is neglected and the adherends are treated as 
membranes, in [3,4] it is assumed that the adherends are membranes and 
the adhesive is a shear spring, in [5-8] the adherends are assumed to be 
plates and the adhesive a tension-shear spring, and in [9-11] one or 
both adherends are treated as an elastic continuum. One should, of 
course, add that by using the finite element method, it is possible to 
treat all three components of the adhesively bonded structure as elastic 
continua. 

In this paper it is assumed that the thickness of the adhesive is 
small compared to the thicknesses of the adherends which, in turn, are 
small compared to the length of the joint. Thus, the problem is formula- 
ted under the following simplifying assumptions: (i) the adherends are 

orthotropic plates for which a transverse shear theory is used, (ii) the 
adhesive is an elastic layer in which the thickness variation of stresses 
is neglected, and (iii) the bonded structure is in a state of plane 
strain, i.e., e z = 0 for the entire structure (Figure 1). The main pur- 
pose of the paper is to show that under the stated assumptions (a) the 
adhesively bonded joint problems can be solved in closed-form and, 

(b) by comparing the solution with that of the finite element method, 
the analytical results thus found are quite realistic. In formulating 
the adhesive a slight improvement is made over the standard tension-shear 
spring model used in [5-8] by taking into account the effect of the 
average in-plane strain e (Figure 1). It should be pointed out that 

A 
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taking into account any material nonlinearity for the adhesive layer in 
solving the problem appears to be quite difficult. The Appendix shows 
the formulation of the problem by using two different approaches to 
account for the material nonlinearity. In both cases the problem is 
reduced to a system of nonlinear differential equations. The questions 
to be resolved in this regard, however, are (a) whether the lengthy and 
tedious effort necessary to solve the complicated nonlinear problem is 
justified for the type of structural problems under consideration, and 
(b) whether a linear (or linearized) rheological behavior of the adhe- 
sive would not be a more important factor than the nonlinear elastic 
behavior affecting the stress distribution and failure in bonded struc- 
tures [12]. 


2. Formulation of the Problem 

Referring to Figure 1 for notation, the. equilibrium equations for 
the adherends which are assumed to be plates under cylindrical bending 
and normal membrane loading may be expressed as 




- T 



dx 


- a 


fl*L = Q 

dx y lx 


h l +h o 


f2x = q . 

dx g 2x 


h 2 +h o 


(la-f) 


where N^ x , Q ix> , (i=l,2) are the stress and moment resultants and 

t * t and o = o are the shear and the normal stress in the adhesive, 
xy y 

Note that in the problem under consideration N xz = 0 = M xz and, aside 
from the shear resultants acting on the boundaries, no other external 
transverse loads are applied to the aaherends. The stress and moment 
resultants are related to the x,y-components of the displacements u., 
v. and to the rotation g. , (1=1,2) by 

I I A 
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(2a-f) 


du 


1T = C l N lx ’ 


diu 

"dx * C 2 N 2x ' 


dv 


dx 


f + 8 lx = Q lx /B l * 


dv. 


IT + s 2x ” Q 2x /B 2 


d 8 lx 

~dT = D l M lx ’ 


d ^2 x 

TT - D 2 M 2x * 


where 


c i 


1 -v . V • r 

1X ,z . B 4 - I H.G., 


h.E. 

1 IX 


1 6 i ixy * 1 


D. = 


12 »-»1x»iz> 


TvfE7 x 
1 1X 


(3a-c) 


Assuming that y-dependence of the strains e x> e^> and y X y in the 
adhesive is negligible; from simple kinematical considerations we obtain 


£ y * < v f v 2> /h o 
h 


Y xy = (u l ' T B lx ‘ u 2 ‘ TT B 2x )/h o ’ 


'xy 

. *'1 l»1 de 1x . 

e x ' dx 2 dx 


du 2 . h 2 d8 2x wo 

“dT + T ~5T )/Z * 


(4a-c) 


where h Q , h^, and hg are the thicknesses of the adhesive and the adherends 
1 and 2, respectively. Since only the plane strain problem is considered, 
the remaining strains in the structure are zero. If E, v, and G denote 
the elastic constants of the adhesive, its stress strain relations may 
be expressed as 


_ l-v-2v 2 V _ T 

e y ' E(l-v) a * 1-v e x » Y xy G * 


(5a, b) 


Through simple eliminations equations (1), (2), (4) and (5) may 
be reduced to the following system of differential equations for the 
adhesive stresses a y = a(x) and x xy = t(x): 
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(6a ,b) 


6, B.0 + 6,fe 


dx^ n dx 


dx _ 


89(7 + 8 3 dS* + 8 4 3T = 0 * 


where 


r h,(h +h,) h,(h+h 9 ) 

1--f (Ct * C 2 ♦ D, + -H-i- D 2 ) . 


_ G / h l D l h 2°2v 
\ 0 “ o / * 


°2 ' TT v ~2 

0 


R - EH-v) r V / h l D 1 ■ h 2°2a J_ fJL + j_n 

8 1 ■ (l-v-2v 2 ) *■ 2(l-v) ( 2 2 ~ h Q B ] + ’ 


Ed-vJCD^Dg) 

8 2 = h Q (l-v-2v z ) * 


„r Mh+h,) h 2 (h o +h 2 } 

S 3 = - [C 1 ' c 2 + J -r- L D 1 ■ - ° 2 


4 °2 ^ * 


8 4 = “ 2h < !(l-v-2v z ) «W>1 - ( h 2 + V D 2^ * 


(7) 


Eliminating a, (6) may further be reduced to 


d 7 x d 5 t d 3 x + dr _ Q 

W + a 5 dx5 + a 3 dx3 + a l dx 0 


where 


(8) 


a 5 — ct^ 8 1 * a 3 " ^*| 8 1 ^2 " a 2 ® 3 ’ 3 1 ^ 1 8 2 " a 2 8 4 ' ( 8 ) 

At the boundaries of the adherends x = + £, the following conditions 
may be prescribed 



N lx U) = ^U) , N 2x (£) = N 2 U) , N 2x (-£) = N 2 (-£) , 

M lx U) = M ] (£) , M 2x (£) = M 2 (£) , M 2x (-t) = M 2 R), 

Q lx U) = Q-, U) > Q 2 x U) = Q 2 U) . Q 2 x U) = Q 2 (-£) . 

( 10 ) 

where N.. (+.£), Q.j(+£), and M.(+£), (i=l,2), are known constants. Note 
that at x=-£ the- boundary conditions for the adherend 1 need not be pre- 
scribed as they must be such that the gross static equilibrium conditions 
of the composite plate are satisfied. 


3. General Solution 


The differential equation (8) is valid for any composite plate in 
-l < x < l which consists of two layers and is subjected to cylindrical 
bending, and membrane and transverse shear loading given by (10). Hence, 
it may be used to solve any bonded joint problem with two adherends. 
Looking for a solution of the form t(x) = e mx , the characteristic equa- 
tion of the problem is found to be 


7 5 3 

m + a^m + a 3 m + a-jin 


= 0 


OD 


From (11) it is seen that 


m 1 = 0 , 

e 3 + a^e 2 + a^e + a^ = 0 

where m 2 = e. Let e-j, e 2> e 3 be the roots of (13) and let 
♦j = , (j=l *2,3). 


( 12 ) 

(13) 


(14) 


Then, the roots of (11) may be written as 


m^ — 0, m 2 - 4* i » m^ — i* 2 > * 

m 3 = ~<p i , m^ = -<j> 2 , my - -<)> 3 > 
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( 15 ) 



where $-| l 2 is real and $ 2 2 and <j> 3 2 are, in general, complex conjugates. 

The general solution of (8) may now be expressed as follows: 

t(x) = K Q + K 1 sinh cosh ^x + K 3 sinh ^x 

+ K 4 cosh 4 > 2 x + K 5 sinh <|> 3 x + K g cosh <fr 3 x. (16) 

where K ,...,K C are the integration constants. From (6) and (16) it 

0 D 

then follows that 


a(x) = “ ““ [(ci-j^i + 4> -| ^ ) ( cosh $-jX + K 2 sinh $-jX) 

+ ( a .|* 2 + 4 > 2 3 )(l <3 cosh 4 > 2 x + K 4 sinh <j> 2 x) 

+ (o-| 4> 3 + <f> 3 3 ) ( K 5 cosh <J> 3 x + K g sinh 4 > 3 x) ] . 


(17) 


It can be shown that expressed in terms of cr and t alone the boundary 
conditions (10) are equivalent to the following relations: 

X 

t dx = N 2 (-£) - N 2 U), 

-l 


l 


cr dx = Q 2 (-£) - Q 2 U) 


(18) 


(19) 


l 


a x dx = M 2 (£) - M 2 (-£) - £Q 2 (£) - £Q 2 (-£) 


-l 


h +h 0 

CN 2 U) - N 2 (-£)] 


( 20 ) 


h n D 


^ t(£) = [C-jN-j (£) - C 2 N 2 (£) y-^-M^(£) ^ ^ 2 (£)], ( 21 ) 


h 2°2 
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•2 r Mh+h,) M h n + M 

-ftp t(£) = jp {[C-j+Cg + 4 d t + 4 D 2^ t 


U) 


h.,D. h 0 D 0 

- “V'Qi^) '^rh W ’ 

[ f~ ( it + ^ * (h 1 D 1 + h 2 D 2 )]a(£) 

0 I £ 

1 -v-2v 2 d 2 . X. rr , VVM D 

• "E'rp^y sf’ ,,(£) fn^r [C i c 2 4 D i 

h«(h +h 9 ) j D-, D« 

- -M- 2 - w w * i£ M i <*> - ^ M 2<« • 


( 22 ) 


(23) 


VBpV * 4TT^7 (h l°l *¥ 2 >]s»i« 

' - if^f & •«« + 2 e; C(h o +h i )D i - <V h 2 > D 2 ><*> 

+ v rc c 1 VW D . MW Di 4 K) 

+ 2(1 -v) LC T C 2 + 4 D 1 4 U 2 J d)F W 

D, 0 9 

= TT " h~ Q 2^ ’ 

0 0 


(24) 


Technically, the integration constants K ,...,Kg may be determined by 
substituting from (16) and (17) into (18-24), which become 


2 K„ 


2 1C 


2 K« c. i\ * «- l\g 

2£K + si nil + -r- 1 sinh <p 9 l + — sinh <fr~£ 
0 4 >i T 1 <|> 2 2 93 j 


= N 2 (-£) - N 2 (£) , 

p- (a 1 +4* 1 2 ) sinh + K 3 (a-,+<i> 2 2 ) sinh <f> 2 £ 
+ K 5 (a-j+<|) 3 2 ) sinh <j> 3 £] = Q 2 (-£) - Q 2 U) » 


(25) 


(26) 
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' ^ Ck 2<V*i 2)(z cosh sinh 

+ K 4 (a 1 +$ 2 2 )(£ cosh <f 2 £ - sinh 4> 2 £) 

+ K 6 ( ai +« 3 2 )(£ cosh $ 3 -e - j- sinh 4> 3 *)] = M 2 (l) 

r ^ 3 h o +h 2 

" 2 U) - *Q 2 (-£) - [N 2 (£) - N 2 (-£)], (27) 

K l*l cosh 4 ^ + K 2 4 > 1 sinh + K 3 4 > 2 cosh <j> 2 Z + K 4 4 » 2 sinh ul 
+ cosh <f 3 £ + K g 4 » 3 sinh 4> 3 * = jp [C ] N 1 (Z) - CgNgU) 

^lOi ^ 7^7 

, (28) 

K o a l + (°i + *i 2 )(K 1 sinh* 1 ^ + K 2 cosh 4 >1 £) + (c^+^MKgSinh*^ 

+ K 4 cosh<^ 2 £) + (a 1 +<t 3 2 )(K 5 sinh4. 3 £ + Kgcosh^-e) 

= " ^ C D-j Q-j CO + h 2 D 2 Q 2 U)], (29) 

3 2 

L (<*!*!+$•, )(x 1 +x 2 ^ 1 ) - o 2 x 3 4. 1 ](K 1 cosh Sinh 4^0 

3 2 

+ [(a^ 2 +* 2 Kx^x^g ) - a 2 x 3 4 » 2 ](K 3 cosh + K 4 sinh * 2 0 

3 2 

+ [(o 1 4>3 + <f 3 )(x 1 +x 2 * 3 ) - a 2 x 3 ^ 3 ](K 5 cosh $ 3 l + Kg sinh ^l) 

3 K“C D 2 M 2W - W*)] • (30) 
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X 4 K 0 - [ — ^ ( x l*l +x 2* 1 ) ‘ x 4“ x 3+l 3( K T s1nh 4»i^ +K 2 c °s h «t> 1 '&) 


“2 

* 

W , 3, , , . 2- 


- [ 1 (x 1 ^ 2 +X 2 <i) 2 ' 3 )-A 4 -X 3 4. 2 t ](K3Sinh4) 2 £+K 4 cosh^ 2 £) 

- (x 1 4> 3 +X 2< i»33)-x 4 -x 3 4»3 2 3(K 5 s in htf.3^.+K 6CO sh( l »3-£) 


- jJ-Cd^U) - d 2 q 2 U)] 
0 


(31) 


where 


M=^(^ + 4 ) '^ (hlDl + h2D2) ’ 


, _ l-v-2v 2 

x 2 — * 

h.(h +h.) M h « + M 

. - v i-r r j. 10 t n tl 0 ±1 n 1 

X 3 “ 2(1 -v) ^ C 1 ■ C 2 4 D 1 5 D9 -* ! 


x 4 " C( h 0 + h l^ D l " < h o + h 2^°2 ^ ^ 

After determining K 0 »...»Kg from (25-31) the adhesive stresses r and a 
may be obtained from (16) and (17), and the stress and moment resultants 
from (1). 

From equations (7) and (6) it is seen that if the adherends have 
the same thickness and the same material constants, then a 2 = b 3 = B 4 = 0 
and the differential equations uncouple. Hence, the problem becomes 
considerably simpler. The elastic results for this case have been given 
in [12] and are recovered numerically from the present formulation by 
selecting the same material constants for the adherends and a very small 
constant for (h-j — h 2 )/h-| - 
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4. Examples 

A single lap joint and a stiffened plate under various loading con- 
ditions shown in Figure 1 are considered as examples. For the five 
problems shown in Figure 1 the constants which appear in the boundary 
conditions (10) may be expressed as follows: 

(a) Lap joint under tension 

N -,(£) = M-j (£) = Q-, U) = N 2 (-£) = M 2 (-l) = Q 2 (-£) = Q Z U) = 0 


h. h,+h« 

N 2 U) = N q , M 2 U) » N 0 (-f - + -^) . (33) 

(b) Lap joint under bending 

M 2 (£) = M q , (34) 

and all remaining eight constants are zero. 

(c) Lap joint under transverse shear 

Q 2 U) = Q 0 » M 2 U) = IQ 0 , (35) 

and all remaining seven constants are zero. 

(d) Stiffened plate under tension 

N 2 U) = N 0 , N \ z (-l) = N q , (36) 

and all remaining seven constants are zero. 

(e) Stiffened plate under bending 

M 2 (£) = M 0 , M 2 (-£) = M q , (37) 

and all remaining seven constants are zero. 


The following material constants and dimensions are used in the 
examples: 

Adherend 1 (a boron-epoxy laminate) 

E lx = 3.24 x 10 7 psi (2.234 x lO^N/m 2 ), 

E lz = 3.50 x 10 6 psi (2.413 x 10 10 N/m 2 ) , 
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G ]xy = 1.23 x 10 6 psi (8.481 x 10 9 N/m 2 ) , 

v ]x = 0.23, h 1 = 0.03 in (7.62 x 10" 4 m) . 

Adherend 2 (aluminum plate) 

E 2 = 10 7 psi (6.895 x 10 10 N/m 2 ) , 

v 2 = 0.3, h 2 = 0.09 in (0.229 x 10“ 2 m) . 

Adhesive (epoxy) 

E = 4.45 x 10 5 psi (3.068 x 10 9 N/m 2 ) , 

G = 1.65 x 10 5 psi (1.138 x 10 9 N/m 2 ) , 
h Q = 0.004 in (1.016 x 10 -4 m). 

The length of the bonded region 21 is 1 in. for the lap joint and 2 in. 
for the stiffened plate. 

The calculated results are shown in Tables 1-4. For three basic 
loading conditions. Table 1 gives the normalized adhesive stresses t and 
a in a single lap joint. The variation of these stresses with l is 
shown in Tables 2 and 3. Table 4 gives the adhesive stresses in the 
stiffened plate. A limited amount of information is also displayed in Figure 2 
in. order to show the trends in the distribution of t(x) and ct(x). In 
the membrane loading of the lap joint, if N-| x (-.£) = N 2x (£) = N Q , then 
to satisfy the static equilibrium conditions additional bending moments 
and/or transverse shear loads equivalent to a couple N 0 (h 0 +h-|/2+h 2 /2) 
must be applied to the structure. For example, if the structure is 
loaded through pin connections, then Q lx and Q 2x at the pins would be 
nonzero and M^ x and M 2x would be zero. In the example considered, the 
equilibrium is satisfied by applying two equal moments at the ends of 
the structure. Needless to say, the results will be dependent on the 
secondary loads applied to the structure to maintain its static equili- 
brium. 
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5. Discussion and the Finite Element Solution 

From Tables 1-3 and Figure 2 it may be observed that in a lap joint 
the adhesive stresses will be concentrated in the end regions of the 
joint provided the length of the joint 21 is large compared to the 
adherend thicknesses h-j and hg* However, as shown by Tables 2 and 3 
for relatively smaller values of l the stresses in the mid-region of the 
joint may no longer be negligible. By and large, the results found in 
this paper are in agreement with those reported in [6]. The differences 
are mainly due to the approximation resulting from the method used for 
the numerical integration of the differential equations in [6] and 
partly due to the differences in the adhesive models used in [6] and 
in this paper' '. Note that because of the smaller bending stiffness 
of the adherend 1, the peak values of a and t at x= -l are greater 
than those at x=£ (see Figures 1,2 and Tables 1-3). This is the phy- 
sically expected result. 

In the solution given in this paper, the peak values of shear as 
well as that of the normal stress in the adhesive are found to be at 
the end points x = +£. This is inherent in the type of the adhesive 
model used in this type of studies. First, referring to the results 
given in [13] with regard to the stress singularities in bonded wedges 
of two dissimilar elastic materials one may observe that if the adher- 
ends and the adhesive are treated as elastic continua then generally the 
interface stresses a and t would have a power singularity at the end 
points On the other hand, if the adherends are treated as elastic 
continua and the adhesive is assumed to be an uncoupled tension-shear 
spring, then the kernels of the related integral equations would have 
only logarithmic singularities and consequently the stresses ct and t 
would be bounded everywhere, including the ends. Furthermore, the 
examples given in [10] and [11] show that in this case the maximum 
stresses are at the end points. In this sense, the condition that the 

- 

1 'For some joint geometries, significant differences were observed 
between the two sets of calculated results. Within the stated 
assumptions, the solution given in this paper is exact. Hence, 
these differences must be due to the numerical integration. 
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shear stress be zero at the end points of the joint as imposed in some 
investigations is not consistent with the tension-shear spring adhesive 
model (see, for example, [7]). Partly to investigate the overall trend 
of the solution and partly to give some idea about the suitability of 
the plate model for the adherends, the finite element method is used to 
solve a sample problem in which the adherends and the adhesive are 
assumed to be elastic continua. 

The problem considered for the finite element solution is that of 
the stiffened plate shown in Figure 3. For simplicity it is assumed 
that the adherends 1 and 2 are of the same material. Material constants 
and dimensions are shown in the figure. The external loads are assumed 
to be either tension or bending applied at the ends of the plate ( i . e . , 
at x= +1 in.). Because of symmetry only one half of the structure is 
considered. Due to the Saint Venant's principle, since the length of 
the extended portion (0.5 in) of the plate is considerably greater than 
its thickness (0.06 in), the details of the distribution of the applied 
loads at the ends have no effect on the stresses in the stiffener and 
in the plate away from the end regions. Rectangular four-node isopara- 
metric finite elements with incompatible modes[14,15] are used in the 
plane strain solution. The final results given in this paper are 
obtained from the mesh assembled with 544 elements interconnected at 
620 nodal points. Because of its extremely small thickness, only one 
layer of finite elements is used for the adhesive. The calculated 
stresses are those at the midpoint of the elements. 

The results calculated by using the finite element method along 
with those obtained from the plate solution given in this paper are 
shown in Figures 4-7. Figures 4 and 5 show the shear and normal 
stresses in the adhesive for a stiffened plate under tension. The bend- 
ing results are shown in Figures 6 and 7. The finite element results 
obtained in this paper are believed to be reasonably accurate. From 
the figures it is seen that the agreement between the two sets of 
results not only with regard to their trends but also from a quantitative 
viewpoint is surprisingly very good. One may therefore conclude that 
in analyzing the adhesively bonded structures with relatively small 
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adherend thicknesses a plate theory properly accounting for the trans- 
verse shear effects may be used with a certain degree of confidence. 

For the case of two orthotropic adherends under plane strain conditions, 
this paper gives the general closed form solution. 
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APPENDIX 


Formulation of a Bonded Joint with 
Nonlinear Adhesive 

Assuming that the adherends are linearly elastic and may be modeled 
as "plates", equations (1-3) are always valid. Again, if we assume that 
the thickness variation of the strains in the adhesive layer is negli- 
gible, the kinematical relations (4) also remain valid even if the 
adhesive is not a linearly elastic material. However, in this case, the 
stress-strain relations (5) must be replaced by a system of nonlinear 
constitutive equations. This may be done in several ways. One approach 
would be to define a strain energy function in terms of the strain invar- 
iants which would, in turn, give the stress-strain relations necessary 
for the formulation of the problem [16]. Thus, if ag, a-j, a Q are the 
strain invariants, defining a strain energy function (j)^, a.|, a Q ) the 
constitutive relations for the adhesive may be expressed as 

> (i >J = x,y,z) . (A.l) 

1J ae.j 

If the material has the same behavior in tension and in compression, 
in its simplest form <f> must have at least six constants. Let us, there- 
fore, assume that 

<J> = A i a 2 2 + A 2 a l + C l a 2 4 + C 2 a 2^ a l + C 3 a 2 a 0 + C 4 a l 2 * ( A - 2 ) 

where 

a 2 = e x + G y + e z ’ ( A - 3 ) 

a l = £ x £ y + Vz + Vz ' 1 (c xy + £ xz + e yz> * 

a o = V/z " ^ € x E yz + £ y £ xz + £ z £ xy ' £ xy £ xz £ yz> > (A ‘ 5) 

and A-j, A 2 , c-j, Cg, Cg, and c^ are the material constants. For the 
plane strain problem under consideration, we have 
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du . d(3 • 

a = f 2 ( u i» v i> 8 ix , -gp ), (i = 1,2) , (A. 11) 

where f-j and are known nonlinear functions. Thus, equations (A. 10) 

and (A. 11) along with the equilibrium equations (1 a-f) and the stress 

resultant-displacement relations (2 a-f) give 14 equations to determine 

the 14 unknown functions t, <j, u^ , v.., N ix , Q ix , and M- x , (i=l,2). 

Since a = 0 the constant c., does not appear in the formulation of the 
0 ^ 

problem. However, even if it were possible to determine the five remain- 
ing material constants with sufficient accuracy, because of the highly 
nonlinear nature of the functions f-j and f 2 » even the numerical solution 
of the problem becomes extremely complicated. 
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Another approach to the formulation of the problem would be to 
assume that the adhesive has a nonlinear elastic behavior of the type 
suggested by Ramberg and Osgood [17] for the deformation theory of 
plasticity. In such materials, the stress-strain relation for normal 
loading is 

e = 5 + o 5 n , (A. 12) 


where c is the strain normalized with respect to cry/E, a is the stress 
normalized with respect to ay, a is a constant (generally 0.02), n is 
the strain hardening coefficient, E is the initial slope of the stress- 
strain curve, and Oy is a constant (the conventional yield strength in 
plasticity). In the three-dimensional case, one may write (see, for 
example, [18]) 

*1J ■ n+v)S TJ + J^o- kk « 14 + fa 5 e n_1 i,j. <U=X, y.z). 

(A. 13) 

where v is a known material constant and 


ij 


= a 


ij 


1 - 

‘ 3 CT kk 


6 . . , 
IJ 


-o O - 

°e " 2 s ij 


s. . 
ij 


(i ,J=x,y,z) 


(A. 14) 


Noting again that for the problem under consideration 

e z = °» e xz = °» e yz = °» a xz = °* G yz = °’ from ( A * 14 ) we obtain 

Ee x = a x " va y " VCT z + \ a ^ ^ (2a x " a y " a z^’ (A. 15) 

Ee y ~ a y ~ V0 X ” va z + J a a e n ^ ^ 2a y " c x " a z^’ (A. 16) 

0 = a z - va x - vcjy + j a 1 (2a z - - a y ) , (A. 17) 

1 T xy ■ (,+v)o xy + 1“ “e""’ °xy ’ (A ' 18) 


where 


a* • 3(^) 2 + (V + A 2 - X ( 


ay 


ay 


a Y 


la x a y + a x a z + cyj z ) . (A. 19) 
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Using now the kinematical relations (4), from (A.15)-(_A.18) it 
follows that 


f h-|Di 

I «¥lx - -T 1 H lx + c 2 N 2x + -¥■ H 2x> 


* <r x - va y - v<r 2 + | (2<yV’z ) * e 


n-1 


f (v r v 2> ■ \ - w x - v,, z * I < 2 vv o z )o ' 


n-1 


0 = o - vo - vo + (2o -o -o )o ^ ^ , 

z x y 2 v z x y' e * 

h h 

2fi^ (u l " T 8 lx " u 2 ' T 8 2x } = (1+v)a xy + T a e 


3a n-1 


(A. 20) 
(A. 21 ) 
(A. 22) 


o 


xy’ 


(A. 23) 


where o , o = o, o and o v .. = t are the stresses in the adhesive layer. 
Thus, (A.20)-(A.23) along with (1 a-f) and (2 a-f) give a system of 
16 equations to determine o x> o y , o z , o xy , u.., v i , e ix „ N ix , Q ix> and 
M 1x , (i = 1,2). 
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Table 1. Adhesive stresses for a single lap joint, 
l = 0.5 in. 



VO , M o = 0 = Q o 

fyo , f 

1 =°= Q n 
0 0 

V°* N <f 0=M o 

x/£ 

t /( N 0 / 2 £) 

o /{ H 0 Zl ) 

t /( Mq / 4 ^) 

a /( M 0 / 4 £*) 

t /( Q 0 / 2 £) 

o /[ Q 0 / 2 l ) 

- 1.0 

- 13.436 

41.191 

345.70 

- 1328.5 

1 

588.10 

- 0.9 

- 3.690 

- 4.295 

83.789 

137.25 


- 68.272 

- 0.8 

- 0.858 

- 1.713 

12.929 

56.239 

■ DeUi 

- 27.432 

- 0.7 

- 0.141 

- 0.555 

- 1.820 

18.880 

11.944 

- 9.184 

- 0.6 

0.011 

- 0.183 

- 3.190 

6.516 


- 3.163 

- 0.5 

0.028 

- 0.063 

- 2.188 

2.357 


- 1.142 

- 0.4 

0.020 

- 0.023 

- 1.233 

0.897 

11.803 

- 0.434 

- 0.3 

0.011 

- 0.009 

- 0.645 

0.356 

11.520 

- 0.175 

- 0.2 

0.005 

- 0.004 

- 0.335 

0.143 

11.362 

- 0.075 

- 0.1 

0.001 

- 0.002 

- 0.194 

0.051 

11.274 

- 0.038 

0 . 

- 0.002 

- 0.002 

- 0.161 

- 0.001 

11.213 

- 0.029 

0.1 

- 0.007 

- 0.003 

- 0.223 

- 0.051 

11.145 

- 0.038 

0.2 

- 0.017 

- 0.006 

- 0.418 

- 0.132 

11.030 

- 0.071 

0.3 

- 0.036 

- 0.012 

• - 0.865 

- 0.292 

10.800 

- 0.147 

0.4 

- 0.078 

- 0.026 

- 1.835 

- 0.632 

10.316 

- 0.313 

0.5 

- 0.166 

- 0.056 

- 3.908 

- 1.362 

9.285 

- 0.671 

0.6 

- 0.356 

- 0.119 

- 8.327 

- 2.939 

7.090 

- 1.443 

0.7 

- 0.764 

- 0.252 

- 17.713 

- 6.369 

2.417 

- 3.113 

0.8 

- 1.641 

- 0.526 

- 37.572 

- 13.767 

- 7.502 

- 6.691 

0.9 

- 3.538 

- 0.851 

- 79.408 

- 23.586 

- 28.512 

- 11.695 

1.0 

- 7.806 

9.693 

- 170.60 

262.06 

- 74.409 

109.72 

j 
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Table 2. Variation of t/N 0 with l for a single lap joint 
under tension (N 0 f 0, M 0 = 0 = Q 0 ) 




x/l 

0.5 

0.4 

0.3 

0.2 

0.1 

- 1.0 

- 13.436 

- 13.436 

- 13.438 

- 13.469 

- 14.111 

- 0.9 

- 3.690 

- 4.846 

- 6.331 

- 8.247 


- 0.8 

- 0.858 

- 1.572 

- 2.796 

- 4.886 


- 0.7 

- 0.141 

- 0.446 

- 1.170 

- 2.843 


- 0.6 

0.011 

- 0.088 

- 0.451 

- 1.635 

- 5.674 

- 0.5 

0.028 

0.010 

- 0.149 

- 0.939 

- 4.601 

- 0.4 

0.020 

0.027 

- 0.034 

- 0.554 

- 3.809 

- 0.3 

0.011 

0.021 

- 0.000 

- 0.359 

- 3.245 

- 0.2 

0.005 

0.011 

- 0.003 

- 0.283 

- 2.869 

- 0.1 

0.001 

0.001 

- 0.024 

- 0.285 

- 2.647 

0 . 

- 0.002 

- 0.011 

- 0.058 

- 0.345 

- 2.558 

0.1 

- 0.007 

- 0.028 

- 0.108 

- 0.459 

- 2.584 

0.2 

- 0.017 

- 0.055 

- 0.184 

- 0.627 

- 2.714 

0.3 

- 0.036 

- 0.104 

- 0.299 

- 0.862 

- 2.944 

0.4 

- 0.078 

- 0.193 

- 0.479 

- 1.183 

- 3.272 

0.5 

- 0.166 

- 0.356 

- 0.761 

! - 1.620 

- 3.701 

o :6 

- 0.356 

- 0.655 

- 1.206 

- 2.212 

- 4.239 

0.7 

- 0.764 

- 1.208 

- 1.911 

- 3.017 

- 4.898 

0.8 

- 1.641 

- 2.230 

- 3.031 

- 4.117 

- 5.698 

0.9 

- 3.538 

- 4.130 

- 4.825 

- 5.635 

- 6.666 

1.0 

- 7.806 

- 7.806 

- 7.805 

- 7.794 

- 7.846 


Table 3. Variation of a/N 0 with L for a single lap joint 
under tension (N 0 f 0, M 0 = 0 = Q 0 ) 



I L fin) 

x/£ 

075 

0.4 

0.3 

0.2 

0.1 

- 1.0 

41.191 

41.191 

41.187 

41.117 

39.711 

- 0.9 

- 4.295 

- 4.292 

- 2.979 

1.487 

12.572 

- 0.8 

- 1.713 

- 2.634 

- 3.818 

- 4.301 

1.169 

- 0.7 

- 0.555 

- 1.093 

- 2.135 

- 3.838 

- 3.180 

- 0.6 

- 0.183 

- 0.444 

- 1.095 

- 2.656 

- 4.458 

- 0.5 

- 0.063 

- 0.183 

- 0.558 

- 1.738 

- 4.473 

- 0.4 

- 0.023 

- 0.078 

- 0.289 

- 1.127 

- 4.030 

- 0.3 

- 0.009 

- 0.035 

- 0.154 

- 0.743 

- 3.477 

- 0.2 

- 0.004 

- 0.017 

- 0.088 

- 0.508 

- 2.956 

- 0.1 

- 0.002 

- 0.010 

- 0.058 

- 0.372 

- 2.515 

0 . 

- 0.002 

- 0.009 

- 0.049 

- 0.302 

- 2.165 

0.1 

- 0.003 

- 0.012 

- 0.054 

- 0.279 

- 1.898 

0.2 

- 0.006 

- 0.020 

- 0.072 

- 0.294 

- 1.702 

0.3 

- 0.012 

- 0.036 

- 0.106 

- 0.343 

- 1.560 

0.4 

- 0.026 

- 0.065 

- 0.163 

- 0.427 

- 1.447 

0.5 

- 0.056 

- 0.119 

- 0.254 

- 0.548 

- 1.322 

0.6 

- 0.119 

- 0.217 

- 0.394 

- 0.702 

- 1.106 

0.7 

- 0.252 

- 0.394 

- 0.605 

- 0.847 

- 0.639 

0.8 

- 0.526 

- 0.689 

- 0.839 

- 0.732 

0.420 

0.9 

- 0.851 

- 0.726 

- 0.283 

0.862 

2.803 

1.0 

9.693 

9.693 

9.691 

| 9.654 

8.075 
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Table 4. Results for a stiffened plate (l = 1 in) 


x/l 

N 0 f 0 , Mf» 

= 0 = Qo 

mmm kri 

EEBEMM 

x/lHM) 

o/iHa/i) 


oMa/W 

0. 

0 

-3.06xl0‘ 7 

0 

-3.01x10"® 

0.05 

-7. 49x1 O' 7 

-3. 99xl0" 7 

-7.36x10"® 

-3.92x10"® 

0.1 

-1 .95x10"® 

-7.33xl0' 7 

-1.92xl0' 4 

-7.20x10"® 

0.15 

-4.34xl0' 6 

-1.51x10'® 

-4. 26x1 0" 4 

-1.48X10" 4 

0.20 

-9.35xl0" 6 

-3.20xl0' 6 

-9.18X10' 4 

-3.1 5xl0" 4 

0.25 

-2. 00x1 0" 5 

-6.84xl0~ 6 

-1 . 97x1 0" 3 

-6.72xl0" 4 

0.30 

-4.29xl0" 5 

-1.46x10"® 

-4.21xl0' 3 

-1 .44x1 0" 3 

0.35 

-9.17xl0" 5 

-3 . 1 2x1 0" 5 

-8.99xl0' 3 

-3. 07x1 0" 3 

0.40 

-1.96x1 O’ 4 

-6. 66x1 O'® 

-0.019 

-6.57xl0" 3 

0.45 

-4.20x1 0‘ 4 

-1 .42x1 0" 4 

-0.041 

-0.014 

0.50 

-8. 98x1 O' 4 

-3.03X10' 4 

-0.088 

-0.030 

0.55 

-1.92x1 O' 3 

-6.44xl0" 4 

-0.188 

-0.064 

0.60 

-4.13xl0" 3 

-1 . 36xlO“ 3 

-0.401 

-0.138 

0.65 

-8. 86x1 O' 3 

-2. 88x1 O" 3 

-0.857 

-0.295 

0.70 

-0.019 

-6.02xl0" 3 

-1.831 

-0.633 

0.75 

-0.041 

-0.012 

-3.907 

-1.362 

0.80 

-0.090 

-0.025 

-8.326 

-2.939 

0.85 

-0.197 

-0.048 

-17.713 

-6.369 

0.90 

-0.439 

-0.085 

-37.571 

-13.767 

0.95 

-0.997 

-0.096 

-79.408 

-23.586 

1.0 

-2.346 

1.307 

-170.60 

262.06 
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Figure 2. The normal (o =a) anc * shear ( T X y = 0 stresses in the 

adhesive for a single lap joint under tension, £=0.5 in. 
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Figure 3. 


The dimensions and material constants for the stiffened 
plate which is solved by using the finite element method. 
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Figure 4. Comparison of the shear stresses in the adhesive for a 
stiffened plate under tension obtained from the plate 
theory and from the finite element method. 



Figure 5. Comparison of the normal stresses in the adhesive for a 
stiffened plate under tension obtained from the plate 
theory and from the finite element method. 
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